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ABSTRACT 

One important aspect to recovering a body falling through the atmosphere, is to 
decelerate and stabilize it. This is usually accomplished by means of a parachute. 
The design of the recovery system necessitates a knowledge of the dynamics and 
loads during parachute deployment and inflation. In many cases, a pitch plane 
analysis will provide adequate information. However, if the body is in a general 
tumbling motion, it is necessary to analyze its motion in ♦’hree dimensions. 

This report contains the equations of motion and a computer program for the dy- 
namics of a six degree of freedom body joined to a five degree of freedom body 
by a quasilinear elastic tether. The forebody is assumed to be a completely 
general rigid body with six degrees of freedom; the decelerator is also assumed to 
be rigid, but with only five degrees of freedom (symmetric about its longitudinal 
axis). The tether is represented by a spring and dashpot in parallel, where the 
spring constant is a function of tether elongation. Lagrange's equation is used 
to derive the equations of motion with the Lagrange multiplier technique used to 
express the constraint provided by the tether. A computer program is included 
which provides a time history of the dynamics of both bodies and the tension in 
the tether. 
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NOMENCLATURE 


The following is a list of variables used in the computer program and in the 
derivation of the equations as discussed in this report. A brief descri, .ion 
and associated units are included. 
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$ 


4 


f ’ 


FORTRAN STANDARD 

A a 

AA(6,4) 

AALPDE(8) 

AALPFE(16) 

AALPME(16) 

AALPPE(8) 

AAM(8) 

AAMD(8) 

AAMP(8) 
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METRIC ENGLISH 

DESCRIPTION UNITS UNITS 

Distance along the longitudinal axis of 

the forebody (X^) from the intersection 

of the body axes to the tether- torebody 

confluence point, positive towards the 

nose (m) ft 


Dummy variables usei to express incre- 
mental velocities of the forebody in 
the Runge-Kutta integration 


(m/sec) ft/sec 

or 

rad/sec 


An array of eight variables signifying 
angle-of-attack of the forebody used 

with damping coefficients deg 


An array of sixteen ■ ariables signifying 

angle-of-attack of the forebody used 

with force coefficients deg 

An array of sixteen variables signifying 
angle-of-attack of the forebody used 

with moment coefficients deg 

An array of eight variables signifying 

angles of attack of the decelerator deg 

An array of eight variables signifying 
Mach number of the forebody used with 
force and moment coefficients 


An array of eight variables signifying 
Mach number of the fortoody used with 
damping coefficients 

An array of eight variables signifying 
Mach number of the decelerator 


AERATO 


Suspension Line AE Ratio 
(AERATO - AE/AE Nylon ) 
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METRIC ENGLISH 


FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

UNITS 

AIPHI 


Number of elements in PPHIE array 
( = 2 to 8) 



AIPHID 


Number of elements in PPHIDE array 
( = 2 to 8) 



AJALPD 


Number of elements in AALPDE array 
( “ 8) 



AJALPF 


Number of elements in AALPFE array 
( = 8 or 16) 



AJALPM 


Number of elements in AALPME array 
( = 8 or 16) 



ARAM 


Number of elements in AAM array 
( = 2 to 8) 



AKAMD 


Number of elements in AAMD array 
( = 2 to 8) 



ALPE 

a 

Angle-of-attack of the forebody 


deg 

ALPPE 

a P 

Angle-of-attack of the decelerator 


deg 

AM 


Mach number of the forebody 



AMAX1 


Larger ordinate of two points on the 
longitudinal added mass versus 
D q log log plot 

kg 

slug 

AMAX2 


Smaller ordinate of two points on the 
longitudinal added mass versus 
D q log log plot 

kg 

slug 

AMAY1 


Larger ordinate of two points on the 
lateral added mass versus 
D q log log Dlot 

kg 

slug 
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FORTRAN 

AMAY2 


AMP 

AP 


AX 


AY 
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METRIC ENGLISH 

STANDARD DESCRIPTION UNITS UNITS 

Smaller ordinate of two points on the 
lateral added mass versus 

D q log log plot kg slug 

Mach number of the decelerator 


a 

P 


Distance along the longitudinal axis 

of the decelerator (X ) from the c.g. 

pb 

to the teth' r-decelerator confluence 

point (m) ft 


Exponent of longitudinal added mass 
equation 

(MPAL = RH000 * BX * DS ** AX) 


Exponent of lateral added mass equation 
(MPAS = RH000 * BY * DS ** AY) 


B b 


BB(5 ,4) 


BX 

BY 

C c 


Distance along the lateral axis of the 

forebody (Y^) from the intersection of 

the body axes to the tether- forebody 

confluence point, positive towards the 

left wing (m) ft 


Dummy variables used to express incre- 
mental velocities of the decelerator in 
tht Runge-Kutta integration 


(m/sec) ft/sec 
or 

rad/sec 


Coefficient of longitudinal added mass 
equation (MPAL = RH000 * BX * DS ** AX) 


Coefficient of lateral added mass 
equation (MPAS = RH000 * BY * DS ** AY) 


Distance along the vertical axis of the 
forebody (Z ) from the intersection of 

D 

the body axes to the tether- forebody 

confluence point, positive up (m) ft 
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FORTRAN 


STANDARD 


CC(3,3) 


CAA(8,16,8) 


CCAP(8,8) 


CCHIP 


CCLL(8,16,8) 


CCLLP(8,8,8) 


CCLM(8, 16,8) 


DESCRIPTION 1 

Axial force coefficient of forebody 

Axial force coefficient of decelerator 

Elements of transformation matrix from 
inertial coordinates to body coordinates 
of the forebody 

A three dimensional array of variables 
signifying axial force coefficients 
of the forebody corresponding to 
AAM(8) , AALPFE(16) , and PPHIE(8) 

A two dimensional array of 64 variables 
signifying axial force coefficients of the 
decelerator with respect to angle of 
attack corresponding to AAMP(l) thru 
AAMP (8) 

cos(X ) 

cos( X p ) 

A three dimensional array of variables 
signifying rol" : ng moment coefficients 
of the forebody corresponding to 
AAM(8) , AALPME (16) , and PPHIE(8) 

A three dimensional array of variables 
signifying roll damping coefficients 
of the forebody corresponding to 
AAMD(8) , AALPDE (8) , and PPHIDE(8) 

A three dimensional array of variables 
signifying pitching moment coefficients 
of the forebody corresponding to 
AAM(8) , AALPME (16) , and PPHIE(8) 


METRIC ENGLISH 
UNITS UNITS 
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FORTRAN STANDARD 
CCLMQ(8,8,8) 


CCLN(8,16,8) 


CCLNR(8,8,8) 
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GER- 16047 
METRIC ENGLISH 

DESCRIPTION UNITS UNITS 

A three dimensional array of variables 
signifying pitch damping coefficients 
of the forebody corresponding to 
AAMD(8) , AALPDE (8) , and PPHIDE(8) 

A three dimensional array of variables 
signifying yawing moment coefficients 
of the forebody corresponding to 
AAM(8) , AALPME (16) , and PPHIE(8) 

A three dimensional array of variables 
signifying yaw damping coefficients of 
the forebody corresponding to 
AAM(8), AALPDE (8), and PPHIDE(8) 


CCN(8,16,8) 


Cv>L (8,8) 


CCNP (8 , 8) 


CCP(3,3) C pij 


A three dimensional array of variables 
signifying normal force coefficients of 
the forebody corresponding to AAM(8) , 
AALPFE(16) , PPHIE(8) 

A two dimensional array of 64 variables 
signifying the pitching moment coefficients 
of the decelerator with respect r.o angle 
of attack corresponding to AAMP(l) thru 
AAMP (8) 

A two dimensional array of 64 variables 
signifying the normal force coefficient 
of the decelerator with respect to angle- 
of attack corresponding to AAMP(l) thru 
AAMP (8) 

Elements of transformation matrix from 
inertial coordinates to body coordinates 
of the decelerator 


* < 
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FORTRAN STANDARD 
CCRIT 

CCY(8, 16,8) 


CDAP 

c.g. 

CGAM C ^ 

CGAMP C v 

*P 

CHIE J 

CHIPE *X 

CLL 
CLLP 
CLM 
CLMQ 
CLN 
CLNR 
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DESCRIPTION 

c/cr = damping ratio s 0.06 

A three dimensional array of variables 
signifying side force coefficients of 
the forebody corresponding to AAM(8) , 
AALPFE(16) , and PPHIE(8) 

Drag area of decelerator 
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METRIC ENGLISH 
UNITS UNITS 


2 2 
(in) ft* 


Center of gravity 

Cos("6 ) 

Cos( ^ ) 

P 

Flight path angle of forebody in 
horizontal plane, measured from X 

axis toward Y axis deg 

Flight path angle of decelerator in 
horizontal plane, measured from X 

axis toward Y axis leg 

Rolling moment coefficient of the 
forebody 

Rolling damping coefficient of the 
forebody 

Pitching moment coefficient of the 
forebody 

Pitch damping coefficient of the 
fore body 

Yawing moment coefficient of the 
forebody 

Yaw damping coefficient of the forebody 


n 
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FORTRAN 

STANDARD 

DESCRIPTION 

METRIC 

UNITS 

ENGLISH 
UNITS 

CN 

C N 

Normal force coefficient of the forebody 


CNP 

C Np 

Normal force coefficient of the decelerator 


COM(20) 


Input variable used to define computer 
simulation - up to eighty figures 



CPHI 


Cos (<S) 



CPHII 

c4. 

1 

Cos 



CPHIPI 

U . 

pi 

Cos ) 

Pi 



CPSI 

c f 

Cos ( f) 



CPSIP 

CS 

C V P 

Cos ( l^ p ) 

Damping coefficient of tether 

N-sec 

lb^-sec 

S 

m 

ft 

CSIGP 


Cosine of one half the apex angle of 
the cone formed by the suspension lines 



CTHE 

C9 

Co (0) 



CTHEP 

CO 

P 

Cos (0 ) 
P 



CY 

s 

Side force coefficient of forebody 



D 

d 

Aerodynamic reference length of 
forebody 

(m) 

ft 

I« (6,6) 


A two dimensional array of variables 
signifying the coefficients of the 
second derivatives in the equations of 
motions 

(kg) 

(kg-m 2 ) 

slug 
or 2 
slug- ft 

DDP(3,3) 


A two dimensional array of variables 
signifying the coefficients of the 
second derivatives in the equations of 
motion of the decelerator 

(kg) 

(kg-m 2 ) 

slugs 
or 2 
slug- ft 

DELSX 


Total suspension line defl ..tion array 

m 

ft 


n 


xi- 
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FORTRAN 

STANDARD 

DESCRIPTION 

METRIC 

UNITS 

ENGLISH 

UNITS 

DELTX 


Total tether line deflection array 

m 

ft 

DLTO 


Initial elongation of tsthev beyond 
unstretched length. DLTO is negative 
if the forebody and decelerator con- 
fluence points are closer together 
than LTO 

(m) 

ft 

DLTX 


Tether deflection component in array 
element (DLX(I) associated with load 
PX(I) 

m 

ft 

DLX 


Effective spring deflection array 

m 

ft 

DMD 


Rate of change of longitudinal added 
mass 

kg/sec 

slug/sec 

DP 

d 

P 

Aerodynamic reference length of 
decelerator 

(m) 

ft 

DPR 


Degress per radian - 57.2957795 



DS 


Parachute diameter associated with SP 

m 

ft 

DSP 


Parachute projected diameter associated 
with DS 

m 

ft 


DSX1 


DSX2 


DSYl 


Larger abscissa of two points cn the 
longitudinal added mass versus 

D q log log plot m ft 

Smaller abscissa of two points on the 
longitudinal added mass versus 

D q log log plot m ft 

Larger abscissa of two points on the 
lateral added mass versus 

D log log plot m ft 

o 



r* 
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FORTRAN 

STANDARD 

DESCRIPTION 

METRIC 

UNITS 

ENGLISH 

UNITS 

DSY2 


Smaller abscissa of two points on the 





lateral added mass versus 





D q log log plot 

m 

ft 

DT 


Integration increment 


sec 

DTP 


Number of integrations between data 





output 



DTFl 


Number of integrations between data 
output when 1‘ i- TDTC 



DTI 


Integration increment when T V TDTC 


sec 

DTT 


Estimated parachute system period/ 12 

sec 

sec 

DYPR 

q 

Dynamic pressure of forebody 

(N/m 2 ) 

lb f /ft 2 

DYPRP 

% 

Dynamic pressure of decelerator 

(N/m 2 ) 

lb f / ft 2 

EE (6) 


An array signifying the nonhomo- 

(N/m) 

ft-lb f 



geneous terms in the six equations 
of motion of the forebody 

(N) 

or 

lb f 

EPL 


Suspension line strain array 

m/m 

ft/ft 

EPS 


Number used to check for incons' >tent 

-13 

equations in PIVERT Subroutine, 10 




EPSI Number used to check if 0 is approaching 

2n+l 6 

a singular point 0 = 2 — If ^ is 

approaching a singular point, the 
accelerations are kept fixed until this 
region is passed. EPSI = O.OOOOObl 
freezes the accelerations if 0 is within 
0.2° of a singular point. 

EPT Tether line strain array m/m ft/ft 
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FORTRAN 

STANDARD 

DESCRIPTION 

METRIC 

UNITS 

ENGLISH 

UNITS 

ETAI 


Number which controls DT if 0 is near a 





singular point. ETAI = 0.00061 sets 
DT = /5, if 0 is within 2° of a 





singularity 



FF(5) 


An array signifying the accelerations 

(m/sec) 

ft/sec 



of the decelerator 


or 

rad/ sec 

FREQP 


Estimated parachute system frequer .y 

1/sec 

I/sec 

FSULT 


Ultimate design factor of safety for 
parachute 



FX 

F 

X 

Generalized force on forebody in 





X direction 

(N) 

lb f 

FXB 

F xb 

Body force in direction of X^ due to 
aerodynamics 

(N) 

lb f 

FXP 

F 

xp 

Generalized force on decelerator in 
X direction 

(N) 

ib f 

FXPB 

F xpb 

Body force in direction of X , 

Pb 

(N) 

lb f 

FY 

F 

y 

Generalized force on forebody in 





Y direction 

(N) 

lb f 

FYB 

F yb 

Body force in direction of Y^ due 





to aerodynamics 

(N) 

lb f 

FYP 

F yp 

Generalized force on decelerator in 
Y direction 

(N) 

lb f 

FYPB 

F ypb 

Body force in direction of Y ^ 

(N) 

lb f 

FZ 

F 

z 

Generalized force on forebody in 





Z direction 

(N) 

lb f 

FZB 

F xb 

Body force in direction of Z^ due to 
ae rodynamics 

(N) 

lb f 
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METRIC 

ENGLISH 

I 

FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

UNITS 

** 

FZP 

F 

Generalized force on decelerator in 





ZP 

Z direction 

(N) 

cr 

>-* 

> 

FZPB 

F . 
zpb 

Body force in direction of Z . 

pb 

(N) 

lb f 


G 

g 

Acceleration of gravity at Z 

(m/sec“) 

2 

ft/sec 


GAME 


Flight path angle of forebody in 
vertical plane 


deg 


GAMPE 

D 

Flight path angle of decelerator in 





r 

vertical plane 


deg 


GLOAD 


Limit design load factor of forebody 




GO 


Acceleration of gravity at earth's 
surface, 32.17 

(m/sec 2 ) 

2 

ft/sec 


HHH 


Altitude below which simulation is 






ended 

(m) 

ft 


ICXO 


Canopy roll moment of inertia about 






its C.M. 

, 2 
kg-m 

2 

slug- ft 


ICYO 


Canopy pitch moment of inertia about 






its C.M. 

. 2 
kg-ra 

2 

slug-ft 


IERSW 


Variable signifying whether or not 
the equations being solved in sub- 
routine PIVERr are consistent 



; 

ILXO 


Parachute lines roll moment of inertia 

2 

. 2 




about its C.M. 

kg-m 

slug- ft 


ILYO 


Parachute lines pitch moment of inertia 

2 

2 




about its C.M. 

kg-m 

slug-ft 

; 

IXB 

\b 

Moment of inertia about axis 

(kg-m 2 ) 

2 

slug-ft 

? 

< 
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METRIC 

ENGLISH 


FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

UNITS 


IXPB 

I . 

xpb 

Apparent moment of inertia about 

2 





X , axis 
pb 

(kg-m ) 

slug- ft 


IXYB 

i . 

Product of inertia associated with 





xyb 







X^ and Y^ axes 

(kg-m ) 

slug-ft 


IXZB 

*xzb 

Product of inertia associated with 






X^ and Z^ axes 

(kg-m ) 

slug- ft 


IYB 

V 

Moment of inertia about Y axis 

b 

(kg-m 2 ) 

slug-ft 


IYPB 

r ypb 

Apparent moment of inertia about 






Y , axis 
pb 

(kg-m ) 

slug- ft' 


IYZB 

I . 

Product of inertia associated with 





yzb 


2 





Y, and Z, axes 
d b 

(kg-.n ) 

slug-ft' 


IZB 

X zb 

Moment of inertia about Z, axis 

b 

(kg-m 2 ) 

slug-ft' 

r 

KS 

K 

s 

Tether spring constant 

(N/m) 

>—• 

cr 

M) 

n 


LS 


Suspension line length 

m 

ft 


LSCL 


Distance along parachute centerline 
between the confluence point and the 
projected diameter plane 

m 

ft 


LI 


Tether length - distance between 
confluence points 

(m) 

ft 


LTD 

\ 

Time rate of change of tether length 

(m/sec) 

ft/sec 


LTO 

hzo 

Unstretched tether length 

(m) 

ft 

~o 

— 05 

M 

m 

Mass of forebody 

(kg) 

slugs 


MP 

m 

Real mass of decelerator 

(kg) 

slugs 

in ui 


p 


o •• 

MPAL 


Added mass of the decelerator along X . 

pb 






axis 

(kg) 

slug? 

— 






r'. 
1 * 
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METRIC 

ENGLISH 


FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

UNITS 

•4 

£ 

MPAS 


Added mass of the decelerator alone Y . 

pb 

or Z , axis 
pb 

(kg) 

slugs 


MPL 

m 

P* 

Apparent longitudinal (X ^) mass of 

P 



- 


decelerator 

(kg) 

slugs 

* 

MPS 

m 

ps 

Apparent side (Y , or Z ) mass of 

pb pb 



- 



dece lerator 

(kg) 

slugs 


NS 


Number of parachute suspension lines 



* 

NT 


Number of tether lines 



-1 

OMETRC 


Option variable: if OMETRC = L, Input 




and Output are in the metric system. 

If OMETRC = 0.0 Input and Output ar-j in 
the English system. 


OMXBE 

^xb 

Angular velocity about X^ axis 

deg/sec 

OMYBE 

Sb 

Angular velocity about axis 

deg/sec 

OMZBE 

" ^ zb 

Angular velocity about axis 

deg/sec 

OPAM 


Option variable: if OPAM = 1., added 



mass of the decelerator f 0; if 


OPAM = 0., added mass of decelerator 
= 0 


OPDA 

Option variable: if OPDA = 1., damping 
moment coefficients of the forebody are 
read in as arrays; if OPDA = 0., damping 
moment coefficients are read in as 


constants 

OPOS 

Option variable: if OPOS = 1., at least 
one of the c.g. offsets or products of 
inertia of the forebody t 0.; if OPOS * 0., 
all c.g. offsets and products of inertia = 0. 
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STANDARD 


DESCRIPTION U1 

Option variable: if OPPLOT = 1., a plot 
tape can be made; if OPPLOT =0., no plot 
tape is made. 

Option variable: if OPPRIN = 1., all aero- 
dynamic coefficient arrays are printed out; 
if OPPRIN = Cl, no aerodynamic coefficient 
arrays are printed out 

Option for automatic DT determination 
(OPDT * 1) 

Option for automatic parachute area 
calculations (OPSP * lj 

Option variable: if OPSYM = 1., the fore- 
body is aerodynamically symmetric such that 

C =C -C =0; if OPSYM = 0, the forebody is 
y m 

not symmetric 

Parachute overinf iation at reefed 
stage (I), (percent/100) 

Parachute overinflation at reefed 
stage 1. (percent/100) 

Parachute overinflat ion at reefed 
stage 2. (percent/100) 

Parachute overinflation at reefed 
stage 3. (percent/ 100) 

Aerodynamic roll angle of forebody, 

0 £ PHIAE n< 180° 

Angular acceleration about axis 
Angular velocity about X^ axis 


METRIC ENGLISH 
UNITS UNITS 


deg/sec 

deg/sec^ 
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FORTRAN STANDARD 


PHIIE 


PHIPI 0 . 

P 1 


POROS 


PPHIDE(8) 


PPHIE(8) 


PSIDDE 


PSIDE 


PSIPDE 


PSIPE 


PSPDDE 


PULAN 
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DESCRIPTION 


METRIC ENGLISH 
UNITS UNITS 


Euler angle rotation about axis 

Aerodynamic roll angle of fore body, 

-180° 0. 4 180° 

Aerodynamic roll angle of decelerator, 

-180° -£0.4 180° 

P 1 

Parachute porosity. Use POROS = 0.15 

An array of eight variables signifying 
forebody roll angle used with damping 
coefficients 

An array of eight variables signifying 
forebody roll angle used with force 
and moment coefficients 

Angular acceleration of forebody 
about -Z axis 

Angular velocity of forebody about 
-Z axis 

Euler angle rotation of forebody about 
-Z axis 

Angula r velocity of decelerator about 
-Z axis 

Angular rotation of decelerator about 
-Z axis 

Angular acceleration of decelerator about 
-Z axis 

Angle between tether and forebody 
centerline 


deg/sec 


deg/sec 


deg/sec 


deg/sec 


Suspension line load array 


O 
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1 

■A * 

FORTRAN 

STANDARD 

DESCRIPTION 

METRIC 

UNITS 

ENGLISH 

UNITS 

V 

PSX 


Total suspension line load array 

N 

lb f 

* 

PT 


Tether line load array 

N 

lb f 

: 

PTX 


Total tether line array 

N 

lb f 

t 

PX 


Effective spring load array 

N 

lb f 

t 

QKAXPB 


Parachute load due to rate of change 
of mass of the parachute times the 
relative velocity, XPBDR 

N 

lb f 


QPHI 

Q « 

Generalized force about X^ axis 

(m-N) 

ft-lb f 


QPSI 

Qt 

Generalized force of forebody about 






-Z axis 

(m-N) 

ft-lb f 


QPSIP 

V 

Generalized for of decelerator about 





-Z axis 

(m/N) 

ft/lb f 

: | 

V - 

QTHE 


Generalized force of .orebody about 
negative line of modes 

(m/N) 

ft/lb f 

y 

QTHEP 

V 

Generalized force of decelerator about 





negative line of modes 

(m/N) 

ft/lb f 


RATIO 


Nondimensional length used in the 
decelerator's body torque expressions 



< 

RE 


Radius of earth - 20,926,435. 

(m) 

ft 

& 

RHO 


Density of atmosphere at Z (1962 
Standard) 

(kg/m 2 ) 

slug/ ft 2 

•~o 

RHOOO 


Air density ratio (RHO/RHOO) 



— oo 
k r-<o 

m uj 

S 

s 

Aerodynamic reference area of forebody 

(m 2 ) 

ft 2 

o •• 

— Lu 
< « UJ 

SCHI 


sin (X ) 




SC1IIP 


sin ( X p ) 



— 

SGAM 


sin Of ) 




o 
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METRIC ENGLISH 


FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

UNITS 

SGAMF 


sin (tf p ) 



SP 


Parachute drag area 

2 

m 

ft 2 

SPD 


Time rate of change of parachute 
drag area 

m^/sec 

ft 2 /sec 

SPHI 

S0 

sin (0) 



SPHII 

S *i 

sin (0^ 



SPHIPI 

^Pi 

sin (0 .) 
Pi 



SPRQ 


Initial parachute area 

2 

m 

ft 2 

SPRl 


First reefed stage parachute drag area 

2 

m 

ft 2 

SPR2 


Second reefed stage parachute drag area 

2 

m 

ft 2 

SPR3 


Third reefed stage parachute drag area 

2 

m 

ft 2 

SPRL 


Parachute drag area associated with 
reefed stage (1-1) 

2 

m 

ft 2 

SPRU 


Parachute drag area associated with 
reefed stage (I) 

2 

m 

ft 2 

SPSI 

s f 

sin (f ) 



SPSIP 

Sfp 

sin< ? p ) 



SSP(16) 


An array of sixteen variables signifying 
aerodynamic reference area of the 
decelerator corresponding to TTIP(16) 

(M 2 ) 

ft 2 

STHE 

SO 

sin (0) 



STHEP 

SO 

p 

sin (0 ) 
P 



T 


Flight time 


sec 

TENS 


Tension in tether 

N 

lb f 

THEDDE 

*0* 

Angular acceleration of forebody about 
negative line of nodes 


deg/sec 2 
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FORTRAN 

THEDE 

THEE 

THEPDE 

THEPE 

THPDDE 

TFI 

TINF 

TINT 

TNINY 

TO 
TO IF 

TOTRO 

T0TR1 

T0TR2 


STANDARD 

0 

0 

9 


P 

9 

P 


0 

P 
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METRIC 

UNITS 


DESCRIPTION 

Angular velocity of forebody about 
negative line of nodes 

Euler angle rotation of forebody about 
negative line of nodes 

Angular velocity of decelerator about 
negative line of nodes 

Euler angle rotation of decelerator about 
negative line of nodes 

Angular acceleration of decelerator about 
negative line of nodes 

Parachute inflation time from Stage (I) 
to Stage (I + 1) 

Time when inflated area first 
equals SPRU 

Tine at start of inflation of 
zeeied stage (I) 

Total time spent in region where 

(1-ABS (SIN(THE)) . LT.EPSI) sec 

Initial time sec 


sec 


sec 


sec 


ENGLISH 

UNITS 


deg/sec 

deg 

deg/sec 

deg 

deg/sec^ 

sec 

sec 

sec 

sec 

sec 


One ' alf the time spent in the over- 
inflation of stage I sec sec 

Time at start of inflation of first 

reefed stage sec sec 

Time at start of inflation of 

second reefed stage sec sec 

Time at start of inflation of 

third reefed stage sec sec 


1 


o 
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FORTRAN 

T0TR3 

STANDARD 

DESCRIPTION 

Time at end . third reefed stage 

METRIC 

UNITS 

sec 

ENGLISH 

UNITS 

sec 

TPD 


Total tether load (tension + damping) 

N 

lb . 

TPDRB 


Component of tether load normal to 
forebody centerline 

N 

lb f 

TPDXB 


Tether load component along forebody 
XB axis 

N 

lb f 

TPDYB 


Tether load component along forebody 
YB axis 

N 

lb f 

TPDZ3 


Tether load component along forebody 
ZB axis 

N 

lb f 

TRO 


Time from "TO" to start of inflation 
of first stage 

sec 

sec 

TR1 


Time from "TO" to start of inflation 
of second stage 

sec 

sec 

TR2 


Time from "TO" to start of inflation 
of third stage 

sec 

sec 

TR3 


Time from "TO" to end of third 
stage (TR3 > TTT) 

sec 

sec 

TTIP (16) 


An array of sixteen variables signifying 
inflation time events 


sec 

TXB 

T xb 

Torque about X. axis due to aerodynamics 
b 

(m-N) 

ft-lb f 

TTT 


Flight time at which simulation is ended 


sec 

TYB 

T . 
yo 

Torque about Y^ axis due to aerodynamics 

(m-N) 

ft-lb f 

TYPB 

T ypb 

Torque about X^ axis due to aero- 
dynamic s 

(m-N) 

ft-lb f 

TZB 

T xb 

Torque about Z^ axis due to aerodynamics 

(M-N) 

ft-lb f 


o 
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FORTRAN 

TZPB 

V 

VP 

VS 

WT 

WTC 

WTCM 

WTL 

WTLM 

WTP 

X 


XBAR 

XBD 

X. 

XDD 

XP 

XPBD 

XPBDR 
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STANDARD 

DESCRIPTION 

METRIC 

UNITS 

ENCLISH 

UNITS 

T zpb 

Torque about axis due to aero- 

dynami cs 

(m-N) 

ft-lb f 


Total velocity of forebody 

(m/sec) 

ft/sec 


Total velocity of decelerator 

(m/sec) 

ft/sec 


Speed of sound at Z 

(m/sec) 

ft/sec 


Weight of forebody 

(N) 

lb 


Parachute canopy weight 

kg*m/sec^ 

lb 


Mass of parachute canopy 

kg 

slug 


Parachute suspension lines weight 

kg*n7sec^ 

lb 


Mass of parachute suspension lines 

kg 

slug 


Weight of decelerator 

(N) 

lb 

X 

Down range inertial axis or displacement 



of forebody 

(m) 

ft 

*b 

Longitudinal body axis or displace- 
ment of forebody 

(m) 

ft 

X 

c,g. offset along X^ axis 

(m) 

ft 

*b 

X^ body axis velocity 

(m/sec) 

ft /sec 

• 

X 

Down rai^a velocity of forebody 

(m/sec) 

it/; iec 

• • 

X 

Down range acceleration of forebody 

(m/sec^) 

ft/sec^ 

X 

p 

Down range disolacement of 
decelerator 

(m) 

ft 

• 

x k 
pb 

X , body axis velo>~ 4 ty of decelerator 
pb 

(m/sec) 

ft/sec 


Velocity of air entering or exiting 
the parachute relative to the 
parachute velocity 

m/sec 

ft/sec 
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FOPTRAN 

STANDARD 

DESCRIPTION 

METRIC 

UNITS 

ENGLISH 

UNITS 

XPBDI 


Parachute velocity, XPBD, at T ■ TINT. 





It is used to calculate fill time 

m/sec 

ft/sec 

XPD 

• 

X 

p 

Down range velocity of decelerator 

(m/sec) 

ft/sec 

XPDD 

• • 
X 
p 

Down range acc leration of decelerator 

(m/sec^) 

ft/sec^ 


Y 

Y 

Cross range inertial axis or displace- 





ment of forebody 

(m) 

ft 


Y b 

Lateral body axis or displacement of 





forebody 

(m) 

ft 

YBAR 

Y 

c.g. offset along Y^ axis 

(m) 

ft 

YBD 

• 

Y b 

Y^ body axis velocity 

(m/sec) 

ft/sec 

YD 

• 

Y 

Cross range velocity of forebody 

(m/sec) 

ft/sec 

YDD 

« • 
Y 

Cross range acceleration of forebody 

r% 

(m/sec ) 

2 

ft/sec 

YP 

Y 

p 

Cross range inertial displacement of 




r 

dece lerator 

(m) 

ft 

YPBD 

V 

Y ^ body axis velocity of decelerator 

(m/sec) 

ft/sec 

YPD 

• 

Y 

p 

Cross range velocity of decelerator 

(m/sec) 

ft/aec 

YPDD 

• • 
Y 

n 

Cross range acceleration of decelerator 

(m/roc^) 

2 

ft/sec 

Z 


Vertical inertial axis or displacement 





of forebody 

(m) 

ft 


Z b 

Vertical body axis or displacement of 





forebody 

(m) 

ft 
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METRIC 

ENGLISH 

FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

UNITS 

ZBAR 

Z 

c.g. offset along Z^ axis 

(m) 

ft 

ZBD 

• 

Z b 

Z^ body axis velocity 

(m/sec) 

ft/sec 

ZD 

• 

Z 

Vertical velocity of forebody 

(m/sec) 

ft/sec 

ZDD 

z 

Vertical acceleration of forebody 

(m/sec Z ) 

ft/sec Z 

ZP 

z 

p 

Vertical inertial displacement of 




r 

decelerator 

(m) 

ft 

ZPBD 

> 

Z , body axis velocity of decelerator 
pb 

(m/sec) 

ft/sec 

ZPD 

z 

D 

Vertical velocity of decelerator 

(m/sec) 

ft/sec 

ZPDD 

r 

* • 
z 

p 

Vertical acceleration of decelerator 

(m/sec Z ) 

2 

ft/sec 
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CHAPTER I - INTRODUCTION 

The system to be simulated is two rigid bodies joined by an elastic tether. The 
forebody may have a completely general shape and mass characteristics, and will be 
free to move with six degrees of freedom (three translational, three rotational). 

The decelerator is assumed to be symmetric in shape and mass characteristics about 
its longitudinal Troll) axis, and will be free to move with five degrees of freedom 
(three translational, two rotational). A frictionless swivel is assumed at the 
decelerator-tether confluence point. Thus the roll motions of the forebody will not 
couple with the decelerator. The tether is simulated Dy a spring and dashpot in 
parallel. Damping coefficients for tether lines are difficult to obtain; but spring 
constants for a tether can be found from experimental stress strain curves. Conse- 
quently; the damping coefficient is assumed constant, while the spring constant is 
assumed to be a function of elongation in the computer program, thereby introducing 
a quasilinear spring. 

r 
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CHAPTER II DERIVATION OF EQUATIONS OF MOTION 


SECTION 1 - COORDINATE SYSTEM 


Figure 1 shows the different coordinate syster s used to derive the 

equations of motion. XYZ is an inertial orthogonal coordinate 

system attached to a flat non-rotating earth. XYZ and X y Z are 

1? J? 

orthogonal axes fixed to the forebody and decelerater at "O" and 

"O " respectively. Coordinate systems XYZ and X Y^ Z^ translate 
P P P P 

with the bodies but do not rotate, always remaining parallel to 

corresponding inertial axes. The displacements X, Y, Z, Xp, Yp, 

and Zp, as measured from the origin of Y.Y 2 , are the six 

translational degrees of freedom of the two bodies. The reference 

forebody body axes, longitudinal (X^) , lateral (Y^) , and vertical 

(Z^) , intersect at "O", the origin of the aerodynamics load system 

of the forebody. The reference decelerator body axes, longitudinal 

(Xp^) , lateral (Yp^) , and vertical (Zp^) intersect at "Op", the 

c.g. of the decelerator. The variables X, Y, Z" are the distances 

frora "0" to the c.g. or che foreoody measured positively in the 

direction of the positive body axes X^, Y^, Z^ respectively. For 

orientation purposes , the reader should position himself as a 

pilot in an airplane. In this position, X^ is positive toward the 

nose, Y fa is positive toward the left wing and Z^ is positive up. 

r^ is the vector distance from the intersection of the longitudinal, 

lateral, and vertical axes of the forebody ("0") to the tether 

confluence point of the forebody 

the c.g. of the decelerator ("0 

of the decelerator; r^ lies along Xp^. 


r 2 is the vector distance from 


the c.g. of the decelerator ("Op") to the tether confluence point 
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SECTION 2 - EULER ANGLE TRANSFORMATION 


In order to specify the angular orientation of a body with reference 
to a non-rotating coordinate system (X, Y, Z) , three succesi"' /- 
rotations are made as shown in Figure 2. The first rotation is in 
the direction, -OZ, such that OX and OY are rotated tnroug * an angle 
ip into Oa and ON respectively. The second rotation is in the 
direction, -ON, such that Oa and OZ are rotated through an angle 0 
into OX^ and Ob respectively. The fin~i rotation is about OX^ 
such that ON and Ob are rotated through an angle <{> into OY^ and 
0Z fa respectively. The three angular rotations (ip, 0, <f> ) specify 
the orientation of the body axes (X^ , Y^, Z^) with respect to the 
inertial axes ( XrY'Z ). Again, from a pilots viewpoint, a 
positive ip is a nose to the right yaw? a positive 0 is a nose up pitch; 
and a positive <j> is a right wing down roll. 


The transformation matrix between the body axes and inertial axes 
is now found by considering one rotation at a time and then combining. 
The first rotation is given by: 


Oa 


Cip 

-Sip 

0 

; | 

r ° x ) 

ON [ 

= 

Sip 

Cip 

0 

1 

OY / 

oz) 


0 

0 

1 


. 0Z J 


( 1 ) 


where Sip = simp and Cip = cosip. 
The second rotation is: 


r > 

ox b 


C0 

0 

S0 

ON 


0 

1 

0 

Ob 
> « 


-S0 

0 

C0 



Oa'! 

- 

I”j 


( 2 ) 
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The final rotation is: 



By substituting Equation (1) into (2) and (2) into (3), the 
transformation matrix [C] is formed 



cipco 

-SipC0 

S0 

j'ox) 


-CipS0S<p+SipC<p 

SipS0S<p+CipC<p 

cos <p 

] 0Y ( 

(4) 

“C^S0C4>-SipS4) 

SipS0C(p-CipS<p 

C0C<J) 

[°zj 



Since [C] is a linear orthogonal ( £ C ik = j, k 


i=l 


jk* 


transformation, its inverse is equal to its transpose. 


(1) 


1,2,3) 

Therefore, 


OX 1 

C<pC0 

-CipS0S<p+SipC(p 

-CipS9C<p-SipSc!> 


r > 
ox b 

OY 

= 

-SipC0 

sips0s<p+apc<p 

S\IjS } C<p-C\pS<p 


0Y b 

OZ 

j 


se 

C0S<p 

cec<p 

1 

oz h 


(5) 


For the decelerator, there is no rotation about the longitudinal 
axis. Consequently, the transformation matrix in (4) is simplified 



by lev-ting <p - 

0. The 

result is 

V- 



[V 

1 

f 

Ct C0„ 
P P 

~ S *p C6 p 

S0 

p 

1 f°V) 

~0 

1 

K 


Sip 

p 

C *p 

0 

1 

0Y P 

Y) Ui 
1 

a •• 

OZ w 
i P&J 


-Cip S0 
*P P 

s V e P 

0 

CD 

i 


OZ I 

l PJ 


- u. 


ft 
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The total angular velocity of the forebody is given by: 


w ■ ¥k " 0^ + (7) 

From the inverse of (3) : 

= C*J - S4>k (8) 

b b 

S. = S*J, + C<frJc (9) 

* bo 

From the inverse of (2' 

k = se? b + ce§ 2 ( 10 ) 


Substituting (8), (9 ) , and (10) into (7): 

m = [-i(Sb)+4>]i, +[-i(C0S4>)-6 (C4>) +[-}(C0C*)+0 (S<J>)]k, (11) 

ODD 

The components of angular velocity for the forebody are: 


xb 

= -ij; (S9)+ 4 > 

• • 

(12) 

XI 

>i 

3 

= -ip (cescj>) -e <C4>) 

• • 

(13) 

0) , 
zb 

Likewise 

= -xp (cec<p) +0 (Sep) 

, for the decelerator, the angular velocities are: 

(14) 

u xpb 

» -i (se ) 

p p 

(15) 

w ypb 

n 

i 

CD • 

tJ 

(16) 

w zpb 

= -J (ce ) 
p p 

(17) 
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The kinetic energy of the system is due to the translational and 

rotational velocities of the forebody and the decelerator. The 

forebody is completely general in shape, and products of inertia 

and c.g. offsets will effect the kinetic energy. On the other hand, 

the decelerator is assumed to be symmetric about the longitudinal 

axis and the aerodynamic loads are referenced to the c.g. Therefore, 

all products of inertia and c.g. offsets are zero. The expression 

(2) 

for kinetic energy is: 


| m[X 2 + i 2 + Z 2 ] + §[I xb u) xb + I yb u» yb 2 + I zb w 2b 2 J 


^yzb^yb^zb + I xzb w xb a) zb + I xyb a) xb a, yt^ 


+ m[X b (a>y b Z-Wzb y > + Y b (u) zb X " u xb Z)+ Z b (w xb Y " V° J 


4v [ v 2] + Ks [ v 2 + v 2] 

+ IF^xpb^xpb + I ypb a) ypb + I zpb“zpb ^ 


(18) 


— o 

PS 


—o 

mu 




In Equation (18) , m^ and m^g include directional mass terms due 

to the air enclosed in the canopy, I b , I ypb / an< * I zpb are 
apparent mass moments of inertia. 
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SECTION 4 - POTENTIAL ENERGY 

The potential energy of the system is due to the gravitational 
potential of both bodies and the elastic potential of the tether. 

I 

V = mgtz + X ( S 0 ) + 7(C0S<j>) + Z(C0C<|>)]+ m p gZ p + ^ K g [ 1 ^ - L Tq ] 2 (19) \ 

•9 

L Tq is the unstretched length of the tether; and L T is the 
stretched length of the tether as given by the geometry of the 
system. Referring to Figure 1: 

L T = 1*2 " *ll (20) 

and ?2 are the vectors from the inertial coordinate system 

(XT? ) to the confluence points of the forebody and 
decelerator respectively. For the forebody. 


P 1 - 

4 -f -► 

Xi + Yj + Zk + r 1 

(21) 

-► 

r l = 

ai b + b Db + ck b 

(22) 

a, b, and c are measured along positive body axes X^, Y^, 
and Z^ respectively. Using the coordinate transformation 
matrix (4) : 


F i - 

[X+a(Ci|»C0)+b(-C^S0S<))+Si|(C4>) + c(-C 4 >Sec<t>-S\pS<p)] i 


+ 

[Y+a(-S^C0)+b(Si^S0S4>+CW)+ c(S^S0C^-C^S<j))J j 


+ 

[Z+a (S0)+ b (C0S4>)+ c (C 0C4> ) ] ic 

(23) 



i 
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Similarly for the decelerator: 


I + *„ 3 + K E 


r_ = a i , 
2 P pb 


Substituting (25) into (24) and using matrix Equation (6), 


? 2 ■ VS'Vp 11 1 * tW S *p C0 p )] ? + tV a o (S9 n>J K 


* P P P P 

L t = [L t -L t ] 1/2 


P p p 


L T = { [X p +a p (.C^ p Ce p )-X-a(C^C0)-b(-Cif/S8S(j)+S4»C4i)-c(-c^S0C4>-S4»S<(>) ) 

+ lY +a (-Sij» C0 )-Y-a (-S\pcd) -b (S\t>SeS<p+C\pC<p ) -c (Si|)S0C<J>-Cij>S$) ] 2 

r r r r 

+ [Z p +a p (S0 p )-Z-a (S0)-b (C0S<j>) -c (C0C0) ] 2 } 1/2 (; 

Define the variables A, B, and C such that: 

L t = {[A] 2 + [B] 2 + [C] 2 } 1/2 C 

Further on in the derivation it will be necessary to know 

the total time derivative of L T and the partial derivatives 
of A, B, and C with respect to the generalized coordinates. 


l t = [AA+BB + C c]/l t 
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£ ■ X +a Tit (-i 


x p +a p % { " s ^p C0 p )+0 p (“C^ p se p ) ]-x-a[^ (-s^c6)+e (-c^se) ] 

• • . 

-b (s^ses<j>-«-c^c 4 >) +e (-c^c 0 s<i>) +<p (-opsec<t>-sips<p)] 

-c [iji (S^S0C<j>-C^S(^f0 (-Ci|>C9C<J>)+$ (C\pSdS<f>-S\pC$) ] 


5 “ Y p + «pI , i'p(“ c, l , p C9 p) + 0p(S^ p S0 p )]-Y-a[i(-CipC0)+0 (S4»S0)] 
-b[i(Ci/)S0S<j)-S^C<j)) + 0 (St|»C0S<|»)+} (S^S0C4>-C\pS<J)) J 

-c (c\psdc<p+s\ps<i>) +d ( sipcect) +i (-siiisest-apcf, j 

£ - Z p +a p [0 p (C0 p ) ]-Z-a [0 (C0) ]-b [0 (-S0S(j>)+5> (C9C4> ) J 
-c[0 (-S0C0i+i(-C0S*)J 


as = 

al 

a? 


lx “ 

"Sy 88 

u a 

Iz 

-l 

as 

3B 

ac 


lx “ 

17 = 

Iz 88 

l 

p 

P 

p 


3A 

3A 

3A 

3A 

*Y “ 

Iz 88 


Iz 

P 

3B 

3B 

dB 

3B 

17 “ 

Iz “ 

’V 


ac 

ac 

ac 

ac 

lx 88 

17 = 

lx 88 
p 

17 

p 
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3A 
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3A 


3A 


3B 

Tip 


3B 

36 


3B 
3 <p 


3B 


3B 


3C 

Tip 


a ( S <iC 6 ) -b ( S ipS 9 S <j>+C ipC 4> ) - c ( S tyS 9 C <J>- C\ps <p ) 

(37) 

a (Ci|>S0)+b (C^C0S<J>)+c (Ci£C9C$) 

(38) 

b(C’;'S9c<})+ sips<p)-c (c\psds<t>-s\pc<ti) 

(39) 

V“ s w 

(40) 

a p ( " c V e p ) 

(41) 

a (Cipce ) -b (CipSQSQ-SipCQ) -c (Ci^S 8C<|)+S^S<f> ) 

(42) 

a {-sijjse ) -b (s^ces<#>) -c cs^cec4>) 

(43) 

-b(s\psec<p-c\ps<p)+c(s\pses<p+apc<t;) 

(44) 

a P ( - c W 

(45) 

a p' s VV 

(46) 

0 

(47) 


3C 

77 


a(-C9)+b(S9S4>)+c(S9C$) 


(48) 


i i 
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7? “ 

b(-cec<j>)+c(c0s.i) 

(49) 

3C 

" 

0 

(50) 

3C 

- 

a p we p ) 

(51) 


SECTION 5 - RALEIGH'S DISSIPATION FUNCTION 


If the viscous damping force is proportional to the velocity 
of the particle at which the force acts, an expression analogous 
to the potential energy of a spring may be used. This function. 


F, is known as Rayleigh's dissipation function, and is defined as 

2 (52) 


( 1 ) 


F = 


1 

I 


n 

Z 

i=l 


c iV 


For this problem, Rayleigh's dampi- ' ‘ < _Ui3idered oniy in 
the tether. 


F 




SECTION 6 - LAGRANGE'S EQUATION 


The Lagrange equation for a non-conservative (aerodynamic forces) 
system with a holonomic (can be expressed as an algebraic 
expression), scleronomous (independent of time) constraint and 
Rayleigh's dissipation function (damping in the elastic tether) 
can be written as: ^ 
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In Equation (54), the term A expresses the generalised 

1.1. 

force exerted by the tether on the i degree of freedom. 
The constraint equation is: 


, _ , - y 1/2 

g « {[A] 2 + (B] 2 + [C]*} - 1^ * 0 


(54) 


(55) 


Q i 


is the generalized force due to the aerodynamics. 



is the force due to damping in the tether. 


Tne Lagrangian is equal to the total kinetic energy of the 
system (Equation (18) ) minus the total potential energy of 
the system (Equation (19) ) . With substitutions from Equations 
(4), (6) and (12) to (17), the Lagrangian is: 


L - 7n.[i 2 + i 2 « 2 ]4m pl IX p (C #p C9 p )+i p (-S* p Ce p )« p (S6 p )] 2 

+ 5*» ps < i* p < s v« p (C V 1 2+ [i P ( - c V V + * P (s W”p (c9 p ) )2) 

+ i b [-J(ces$)-9 (c4>)) 2 +j i zb r-i(c0c*)+e(S6)] 


• • 


” {I yzb (C0S<>) _e (c<|,) ] ( cec< fr) +e ( s *> 1 +I x2 b [ “* ( se > + ^l [-* (C6C*)+e (S*) ] 


+ I xyb (-*(S0)+4>] [-*(C0S*)-0(C*)] } 
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+ 5 *xpb >-»p 2 * 1 J ypb [ ®p 2 + *p 2 «* » p > I 

+mt [x(Cij)Ce)+Y (-s^ce)+z (se) ] [z (-^(ces^j-e (c<i>) )-Y(-i-(cec<j»)+e (s*) )] 

+ [x(-oi>ses<M-sij)C4>)+Y (s<ps0s<j>+oj>c<j>)+z (ces<j>) 3 [x(-i(cec<i»)+8 (s$ ) )-z (— i (se)+^) ] 

+ [X(-CiJ>S 0 C<J>-StJ/S 0 )+Y (SiJ/S 0 C$“Cij>S$)+Z (C 6 C 4 > ) ] [Y (~i (S 6 )+J>-X(-J (C 0 S<fr )-6 (C<J>) )] } 
**mg[Z+X(S 8 )+Y(C 0 S<t)+Z(C 0 C<|i)]-m p gZ p - 2 K s [ 1 ^- 1 ^]^ ( 56 ) 

Note: I ypj 3 = I 2 pb ^ ue to decelerator symmetry. 

SECTION 7 - GENERAL EQUATIONS OF MOTION 

Equation ( 56 ) displays all of the generalized coordinates 
explicitly except those appearing n L^,. The terms to be 
substituted into Equation ( 54 ) are now developed. 

X Equation 

ir (^r) * X{m}+^{m(X(-SipC0)+Y(S<|;S0S(()+Ci|/C(j>)+Z‘(Si|>S0C«|>-Cij)S«|;)3 } 

QZ 3 X 

+ 0 {m [X(-Cij<S 0 )+Y(-C^C 0 S«}i)+Z (-Ci^C 0 C 4 >) ] } 

+V{mlY(-OJ/S 0 C<f>-Si|>S$)+Z (Ci^S 0 S<|>-S^S<j>) j } 


+<!> 0{m[X (2S>J»S8)+Y (2SipC0S$) + Z (2Sij;C0C0) ] } 
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+ (m [7 ( 2S 4)S 0 C <fr- 2C\(/S $ ) +Z ( - 2S <^S 9 S <j>- 2Cij)C* ) ] } 

+e}{m[y (-2C(j)cec4>)+?'(20j)ces<<)) ] } 
+J» 2 {mtX(-C^C0)+7(C^S0S<j)-S^C^)+Z‘(C^SeC<J)+S^S(|)) J } 
+0 2 {iri[X(-C4»C0)+Y(C^S0S^)+Z(C^S6C(j!)] } 

+<|> 2 {m[Y(C«S0S$-S(JiC<f>)+Z (CiJ>S0C$+S(|iS4>) ] } (57) 

(58) 

(59) 

(60) 


£r-(— ) = Y{m}+V{m[X(-CiJ>Ce)+Y(C(JiS0S0-SipC4>)+Z(C(|iS0C<t>+SiJ)S<|>)]} 

3Y 

+0{m[X(S(j»S6)+Y (S'<pCdS<p) +Y (Si^C0C<f>) ] } 

+V{m [Y (Si^S 0C4>— CipSif* ) +Z" (-S(|)S0S<j>-Ci|>C<J>) ] } 

+ J»0 {m (X (2C^S0 ) +Y (2C\pC6S$) + Z (2CiJ>C0C<t> ) J ) 


w = 0 

3g -A 

n - l t 
4 - o 

3X 

Y Equation 
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+5|><$>{m[Y (2Cij)SeC<|>+2S5|)S(Ji)+2(-2C^S0S<j)+2S^C<|)) ] } 
+eJ{m{Y(2S^C0C«J))+Z(-2St|)C0S^)] } 

+ji 2 {m [X (S(j;C0 ) +Y (-S5j>S0S4)-Cii>C$) +Z (-S(J)S0C<M-CiJ/S<M ] } 

+0 2 {m [X (S(|»C0 ) +Y (-S^SOStj) ) +Z(-S(|)S0C(|)) ] } 


+<|> {m[Y(-S(JiS0S$-CipC<|>)+Z (-S*S0C<|>+C(|/S<f>) } 

(61) 

3L 

w = 

0 

(62) 

ia = 

3Y 

- B 

S 

(63) 

3F 

3Y 

0 

(64) 

Z Equation 


d / 3L i 
at 3Z 

= Z{m}+0{m[X(C0)+Y(-S0S4i)+K-S0C(j>)J } 



+${m[Y(C0C<J>)+Z(-C0S<t) 1 ) 



+ 04>{m[Y(-2S0C^)+Z’(S0S<|)) ] > 



+ 0 2 {m[X(-S0)+Y (-C0S<(>) +Z’(-C0C<J>) ] } 



+i 2 {m tY(-C0S<f»)+Z(-C0C$)] } 

(65) 
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3X P 

+y P {tm P^ m PS lt -5 S2 *P c2e p I 


+ V t W 1 [ T c *p S 20 p J > 


( Vp ( "V»-V"’ 52 *p c V 11 

+ Vp<'v-v [ - c V 29 p !> 


+Y P V ‘VV 1 ‘• c2 *p c V } 


+Y p e p t[ V- m ps H 3 s2 V 2e P n 


+ Vp { IVV 1 ‘*3 S V 29 P ] 1 


p p 


+2 p e P {Im pr , "p 3 HC *p C2e F 1) 


I 
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3L 

7x_ 


(70) 


3£ A_ 

■5x_ L_ 


(71) 


4 - 0 

3X_ 


(72) 


Y Equation 

P 


d ,3L 

at 


(ff-) - ‘y } 

p 

+ V p (m pi tS 2 * p C 2 0 p J + m ps tc 2 i> p +s 2 i|> p s 2 e p : } 


+S p { V“k.. ,, T s V 29 P n 


^VW f-C2*pC 2 8 p J} 


+i p 5 p l |m pr"'ps 1 1 5 S2 *P S26 P :I 1 


p p 


+i p^p {tm pr m P s 1 IS 2 Vp ]} 

+ Vp { ‘W I-S 2 V 2 e p n 
+ V P { ‘w 1 t-^V 2 V 1 
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3L 

W " 

P 

0 

(74) 

■ 

B 

l t 

(75) 

3F 

3Y 

P 

0 

(76) 

Z p Equation 



at - V'v v 1 [ i c *P S 2 e p 1} 


+ v'w 11 ** s v 2e p n 


p p 


+ Z {m 9 [S 2 0 ]+m „ [C 2 e ] } 
P P* P PS 1 P J 


+ i p*p<i"'pr"'ps I ‘-5 s *p S29 p> ! 


+ X 0 { [m 9 -m ] [Ct|» C20 ] } 

P P P* PS J *P P J 


+ Vp"VV i-2 c *p S29 p ]) 

+ Vp'ivV^V^V 1 

+ Vpii-pryiszyi 


(77) 
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(78) 

i£_ = 

3Z 

P 

C_ 

l t 

(79) 

3F 

3Z 

P 

0 

(80) 

y Equation 


d ,3L 

at ( 7f 

3 ip 

) = 

V {l xb (s2a) + [I yb (s2 ^ )+I zb (c2<}>):1 (C 2 9)+I yzb (-C 2 0S2<(>) 



+ [I x 2 b (C ^ )+I xyb (S< * >)] (“ s20 ) } 


+5 ‘ 1[{I yb- I 2 b ) (7 s 2*)- I y 2 b (C 2 * )J(C 9 )+tI x 2 b (S » ) - I xyb (C * )HS0)> 

+*{i xb (-se) + [i x!!b (c*)-fi jtyb (S4 1 )](ce)} 

+x{m[x(-siiic0)+y (s^ses<j)+c^c4i)+z(s^sec<})-'C(|iS4))] } 

+Y{m [-X (Ci^ce ) -Y (-C^SeS $+S(JjC4) ) -Z (-CtpSGC \pS<t>) ] } 

+$!§{ [I xb “I yb (S 2 <P)-I zb (C 2 0)+I yzb (S24>) ] (S20 ) - [I X2b (C4>)+I xyb (S4») ] (2C20) } 
+Ui t,i yb -I zb ’)(S2((»)-I yzb (2C2^)] (C 2 0) + II xzb (S$)-I xyb (C^)] (S20)} 

+6$ { t- I xb -^ yb “ I zb ) (C2<I>)+I y 2 fc (2S2< ^ ) } 
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+ e2{[(l yb -i zb ) (-• jS2(Ji) +Iy zb (C 24 >) ] (S9)+I x2fc (ceS({.)-I xyb (C8C4))]} 

+ ^ 2{l xyb {CeC$) " I xzb (C6s ^ )} 

+xij; {m t-x {C\J)C0 ) -Y (-Ci|>s 0S(}!+SijjC!f> ) -K-C^S0C(j>-S^S(t» ) ] } 
+X§{n[X(S4)S9)+Y(S^CeS<Ji)+Z(S^C0C<|>) ] } 

•i-x* {m [Y (S^S9C4)-C^S4» ) “'ZtS^SOS^+Cii/C^ ) ] } 

+Y<j/{m [-X (~S\p ce ) -Y (SipSdS<p+C4>C<f>) -"Z (S\pSdC<p-C\pS<}>) ] } 
+Y0{m[X(CipS0)+Y(C^C9S4>)+Z (C\pC6C<p) ] } 

+Y$ {in [-Y (-C\pSdC<p-SipS<t)'j +Z (-C^S 0 S 4>+S^C«f) ) ] 

= xij/{m [-X (C^C0 ) -Y (-C4>S0S$+Sy>C4>) -7 (-C\pSdC<p~SipS<p) ] } 

+xe {m [x (S4ise ) +y (s\j»ces*)+z (sij>cec$>) ] } 

+X4> (m [Y (S ipS 0C<J>— CipS cj> ) — z” CSi^S 0 S 4>+C4»Ccj> ) ] } 

+Y${m [-X (-SiJjC 0 ) -Y (S^S0S<£+C^C<J>) -Z” (S^S 9C4>— Ci^S4>) 3 } 

+yS {m [x (apse) +7 (cipces ^ ) +1 (apcec<t>) ] 

+$${m[-7(-CiJ>S0C4>-SipS4>) + Z^ (~C\pS3S<p+SipC<p) ] } 

{ * If + 1 If } /l t 
0 



C81) 


(82) 

(83) 

(84) 
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at 'rf 1 ' ♦'{[U b-^b* (jS 2 *)-i y2 b ( «*H( c e)-[r b <e*)-i xlb <s*)](S9)} 

0 U 

+6{ I y b( c2 * ) + I zb (E 2 *)t I yzb (S2* ) } 

‘ilv lc( l-'< I b ISH1 

+X{m [X (~C\pS6) +Y (-C^C0S4> ) +Z (-ClJ/C6C<J>) ] ) 

+Y {m [X (S^S 0 ) +Y (S ipC9 S 4> ) +Z (S \^C0C4> ) ] } 

+*z{m[x(c8)+Y(-ses<M+z(-sec$) ] > 

+Ve{ [ <I yb -I 2b ) (-§S24.)-hI yzb (C2*> ] (S9)+ [I^h (C4>) -t xzb (S<|>) ] (-C0) > 
+ HU(I yb -I zb i (C2$)+I yzb (2S2*)] (C0) + tI xyb (S<J.)+I xzb (C<{.)] (SO) } 

+ e^{[I yb -I zb J (-S 24 >)+I yzb (2C2<j>) } 

+ i 2{ - I xyb <s * ) - 1 **b <c *> ) 

+X^{m[xfe^Se)+Y(S^C0S<())+Z'(SipC0C(|>)] } 

+X§ {m [X C-Ctpco } +7 CCip$es<p ) +7 (C^s 0C<J> } ] } 

+X$ {m [7 (-C^COC* ) +7 Cci^ce S<p ) ] > 


+i {m [7 CC^S 0 ) +7 CCi^C 0 5 41 ) +7 (C4/C0C<f> ) J } 


1 
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+y& {m [x (s;ij/ce ) +7 (-siiises^ ) +7 (-si|>s0c<j» ) ] } 

+Y5»{miY(£4»cec4))+z (-s4>ces<t>) ] ) 

+Z§ {m [-X (S0) -Y (C0S$) -Z (C0C<j>) ] } 

+Z$>{m[Y (-S0C<t) + Z (S0S$) 3 } (85) 

H - WUI zb -I yb ](|seS2*) + [-I xyb (C*) + I x2b (S*n(C0) + !I yzb ](S0C2*)} 

+ni[i xh ] (-c 0 )+[i x2b (c<j>)+i xyb (s<j>)] (-S0)> 

♦* a{ ‘ I xb- 1 y b ,sa «- 1 «b <ca *> +I yzb (S2({))] (^S20)-[I xzb (C<t>)+I xyb (Sd))) (C20) } 
+X4<{m [X (S^S0) +Y (Si|»C0S$)+7(S^C0C*) ) } 

+X§ {m [X (-C^C0) +7 (C^S0S(J))+Y (C^S8C$) ] } 

+X$ {m [7 (-Cil>C6C<p) +Z (CxpC0S<J> ) ] } 

+Y4»{m [X(C<J>S0)+7 (C\pCQS<t>) +7 (C^C9C<|)> ] } 

+Y0 {m [X (S^C0 ) + 7 (-SipSeSQ ) +7 (-SipS 0C4> ) ] } 

+Y$ {m [7 (S\J>C0C<|> ) +7 {-Si|>C0S<f> ) ] > 

+Z§{m(X(-S0)+7 (-C0S<J>)+7 (— C0Ccf>) ) } 

+Z${m[7(-S0C4>)+Z <S0S4>) ] } 

-mg[X(C9)+Y(-S0S<j>)+Z(-seC<j>)] (86) 


- 24 - 


- ID- 15(7 -71 } 
EF: EOI 380 


) 


J 


U) Ct 


GOODYEAR AEROSPACE 

COMOOTIO* 



If 
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7 ? 3 ^ 
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(87) 


3F 

—■ = 0 

39 88) 


<t> Equation 


d /3L. v,. , 

3t ( a ^ } = * xb ( ' S0)+I xyb (CeS<J,, + I x 2 b (CeC(J>)} 

+ ® {l xyb (C ^ ) " I xzb (S * ) }+ ‘^ I xb } 

+X{m[Y (-CtySdC<p-S4iS<p) +Z (CipS8S<p-S\pC<p) ] } 

+Y{m [Y (StpSdCip- CipS<fi)+Z(-Si}>SdS4-C<pC<t>)] } 
+Z{m[Y(C0C<J>)+Z(-CeS<))) ] } 

+4>e {l xb (-ce)+ [I xyb (S4>)+I xzb (C<j>) ] (-S6) } 
rri xyb (C ^ ) “ I x 2 b (S< ^ )I (C0 ) } 

+5^-i xyb < s +)-i xzb (c»)> 

+X^{m [Y (8i/;S0C<^“C^S^)+2" (~S](/SdS<t>-CipC<p) ] } 
+X0 {m (Y (-C^C6C()>) +Z* (CipCQStp) ] } 

+X${m[Y (C4>S8S<p-Sii/C<p)+Z (CipSdCQ+SipSip) ] } 
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+Y§{m[Y(SipC0C<i>)+Y'-SipC0S<j>) ] } 

+Y^{m[y^S4-S0S$-Ci^C(t>) + Z(-S^S0C4)+C^S^)] } 

+Z0{m[Y(-S0C4>) + Z(S0S<j>) ] }+Z$ {m [Y (-C0Scp) +Z(-C0C<p) ] > (89) 

|| - ie(t(i yb -i zb )(C2*)+i y2b (2s2»)](ce) + ii xyb (s»)+i xzb (c*)](se)) 

+ 54 { t X x yb ^ C ^^“ I X 2 b t C0 ^ } 

+S}t-i xyb (s*)-i x2b (c * ) ' 

4 2 f[(I yb -I zb ! (^2*)-I y2b (C2»)HC 2 0) + [I x2b (S+!-I xyb (C<M] (^S2en 

+s 2 {(I yb- I z b> ( - 7 s 2 * )+I y Z b (c2 * )} 

+X^{m [Y CS^S0C4>-C^S4>) +7 C-SipS 8S<p-OpC<p) ] } 

+X0{m[Y(-CipC0C<p)+7(CipC0Sct>) ] } 

+x5> (m [Y (Ops 0S (j)-S^C(fi ) +Y CC^S 0C<t)+Si^S <|> ) J 

+Y$ {m [Y (CtpS 0C<p+S;pS <p ) +Z (-OpS 0 S 0+SipC<j> ) ] } 

Hi +Y0{m{Y(S^pC0C<p)+Z(-SipC0S4>) J } 

wO 
m u 

Sul +Y4){m[Y(-SipS0S<j)-OpC<p)+Z (-SipS0C(p+OpS4>)] } 
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ib Equation 
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0 Equation 


d , 3L . 

U ^ A A 


1 v.9 

ypb p 


" X p 2 Um p( - m ps J[-lc 2 4,pS20pl) 

+{ p 2{ 'W '-5® 2 V 2 V 1 

«p 2 { 'VV 11 3 S 29 P 11 

+ Vp { ' W 1 tj 62 v 52 V 1 


+ X Z { [m. -m ] [Cip C26 ] } 
p p pi ps p p 

+Y 2 { [m -m l [-S^ C20 l } 
p p pi ps P P 

2 { [I , -I , ] [is29 ] } 
r p xpb ypb 2 p 
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L t Equation 


3t< 

<7X1 | 

• * 

3L r 

= 0 

(101) 

3L 

3l t 

= 

- K s < V L TO> 

(102) 

3l t 

= 

-1 

(103) 

3F 





= 

C s £ T 

(104) 


Substituting Equations (101) to (104) irto (54) yields: 

A = -[Kgd^-^J+Cg^J (105) 


The value of A in Equation (54) is now defined and is expressed 
in '■erms of the eleven generalized coordinator and their time 
der. natives. The eleven simultaneous, nonlinear, 
coupled differential equations of motion are then written as: 

1. X Equation 

x{m}+ii*{m[x c 12 +y c 22 +z c 32 H 

+e{m[X(-C^S0)+Y(-C'4>C0S$)+Z {-C*C9C4>) j }+$ {m [YC^-ZC^J } 

- ^0{-2m[X(S^S9)+Y(S^C0S<|»)-i-Z(SipC0C<{>)] }+^{-2m[YC 32 -ZC 22 ] } 
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+0|{2m[Y(C4<CeC4>)-Z(C^C0S<Ji)] }+i 2 {m[XC 1]L +YC 21 +ZC 31 ] } 
+§ 2 {m [XC n +Y (-C<^S0S<j>) +Z (-C^S0C» ) ] }+$ 2 (m [YC 21 +ZC 31 J } 


+ IK s ( V L TO )+C sV tA/L T J+F x 


(106) 


2. Y Equation 


•V . •• 


Y{m}+4;{-m[XC 11 +YC 21 +ZC 31 ] } 


+ 0 {m [X (S^S0 ) +Y (St^C0S<j>)+Z (S^C0C<J>) ] }+${n [YC^-ZC^J } 


= ^9{-2m(XfC^S0)+Y(C4/C0S4>)+Z(C^C0C(t))) }+^{2m[YC 31 -ZC 21 ] } 

+0$ {-2m [Y (S^CQC^) -Z (S<fiC0S<f>) ] }+^ 2 (m [XC 12 +YC 22 +ZC 32 ] } 


+0 2 {m[XC 12 +Y(S^S0S4»)+Z (S^S0C$) ] }+$ 2 (m [YC 22 +ZC 32 ] } 


2 r_ r 57 /" 


+ [k s (l t' l to )+c sS 1ib/l t] + t 


3. Z Equation 


# 0 

Z{m}+0*{m[X(C0)+Y(-S0S$)+Z(-S0C<)))] }+4>{m[YC 33 -ZC 23 ] } 


= 0${2m[Y(S0C<M-Z <S0S$) ] } + 0 2 (m [XC 13 +YC 23 +ZC 33 1 } 


+$ 2 {m[YC 23 +ZC 33 ] }+[K s (L t -L to )+C s L t ] [C/L T ]- mg+F 2 
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4. X p Equation 


Vv t '- 2 V 2 e p )+m ps ls 2 V c 2 V 2 V } 


+ V VV [-jS2* p C 2e p I 1+Z p { ‘W 1 [jC*pS20 p J ) 


P P P P* PS J 


r P P 


1 < W { Vp V + Vp < c2 *p S2 V + Vp l« VV 

+ Vp I-^V^V+Vp l ¥' V^V+Vp t- c *p c 2 «p] > 


• |, S l V l I 0 1,C sV 1A/ V + P xp 


(109) 


5. Y Equation 
P 

V vv ‘-^v'y ,+ Vv ts V 2 v + v ic V s V 2 v 1 

+ ' z 'p llm P r m ps lt - 5 s V 29 P n 

- vw fVp^y^^Vp'-^^vvp^V'i. 1 


+Y n e [S 2 4 i S29 ]+Z $ [idl* S20 ]+Z 0 lS<Ji C20 ]} 
P P P P P P 2 y p p J p p l *p p J 


“ [K (L_-L rn/> )+C L_] [B/L’+F 
s T TO' s T J 1 ' T- yp 


(110) 
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6. Z p Equation 


V Im p r m ps ] l T- V 2 V )+ V [m pt‘ m ps If '¥*p S29 p I 1 

+ ' z p""pj (s2e p 1+m ps [c2 y * 

- V^. H WH S2 V + Vp [ ‘ C V :29 P 1+ Vp t ^V 2 V 
+ Vp IS V :29 P I *Vp I ' S2 V , 


"^(Lr-LjQJ+CgLj] [C/L T ]-m g + F 


zp 


(111) 


7. ill Equation 

? {1 ^ ( s 2 e) + [iyb ( s 2 * ) + l zb( c %)](c 2 e )+ i bC-c 2 ® 82 #) 


+ [I xzb (C * )+I xyb (S * )] ( ‘ S26) }+0{[(I yb~ I 2 b ) ^yzb (C2 * ] 1 (C0) 

^W^xyb^ 3 {S0) )+ ^- I xb C 13 +I xyb C 23 +I xzb C 33> 
+X{m[XC 12 +YC 22 +?C 32 ] }+y{-m[XC 11 +YC 21 +ZC 31 ] } 

= ^Ul xb -I yb (S 2 4.)-I zb (C 2 ^)+Iy 2b (S2<j))] (-S26)+[I xzb (C*) 

+1 xyb (S * )] <2C20)}+^{[(I yb -I zb ) (S2<fr) -I yzb (2C24>) ] (-C 2 6) 

+ ^xzb (S *> ** I xyb (C *> J ( “ S20) }+0 * { [I xb‘ ^ I yb~ I zb ) ^^yzb (2S2 *> ] (C6) } 
+ '* 2 < t ^yb' 1 ^ 2 ** ^yzb (C2 *> ’ < S9 > 'Wza^xyb^ > 
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{ - I xyb C 33 + ‘xzbSj 1 

-'WW 6 .** 1 ' 5 If + 5 + 0* 


(112) 


8. 8 Equation 

x{m (x (-ojise ) +y (~c4>ces<p ) +z (-ai>cec<p ) j } 

+Y{m[x(sij/se)+Y ( sii/cesfjo+z - (s^cec^) ] } 

+Z{m[X(Ce)+Y(-S0S4>)+Z (-S9C<t>) ] } 

+*( I dyb-Izb^^’-Iyzb (C2 »>’ 1 < C6 > + tl*yb < C *>- Z XXb (S ♦> 1 <' S9 ' > 

+ 0{l yb (C 2 *)+I zb (S 2 *)+I yzb (S24,)}+i'{I xyb (C*)-I x2b (S*)} 

- M<I I x b +<I yb' I zb* (C2l, ’ ,+I yzb (2S2 ' t)! (-C9)+U xlb (C*)+I b ts+)] (-2S6!) 
♦ 5 » {(I yb- I zb> (S2 +>- I y zb (2C2 + )> 

+ * 2{ ^xb'V (s2 *> -^b (c2 *> +I yzb (s2 »> 1 'I 8128 * + l r xzb (c »> 

+I xyb (S<J>)] <" c20 >> + ^ 2{l X yb (S< ^ )+I xzb (C<|,)} " rng IX(C6)+Y(-S0S<|))+Z(-S8C4.)] 

“ t R s * l t“ l to* +c s^ Te + B Te + c + Q e (H3) 


9. <p Equatio n 

X{m[YC 31 -ZC 21 ] }+Y{m[YC 32 -ZC 22 ] }+Z{m [YC 53 -ZC 23 J } 
+ * l - I xb C 13 +I xyb c 23 +I xzb c 33 1+9{I xyb (c *>- I xzb (S ^ )+ * (I xb> 


2 
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* 5{II xb +(I Yb‘ I ^ ) ( C 2 »> +I yzb <2s2 * )1 (c9) + II xyb (S4,)+I xJb (c * H (2S0)} 
+* 2 {[(i yb - I lb )<5 s2 * ) - I yzb (c2 * )I(c29) + II xxb (s * ) - I xyb (c * ) ^ s29)) 

*® 2 { <I yb _I zb ) (-fS 2 *> +I yz b (C2 * ) )-i"q[V C 33-Z c 23 ] 


. , rT 3A x r iB . = 3 C]/L t + Q. 
- I K s (l t- l to )+c s l t 1 ia T? + B 9? +c * 


(114) 


10. Equation 


* {i . (s 2 e )+i „ H (c 2 e_) } 

y p l xpb p ypb P 


■= y p < ^xpb^ypb' <- s29 p> ,+ ‘W 1 < Vp [ - C 2 V 20 P ] 

+y p i-H s 29 p 1 + Vp i -5 c * P s29 P i+i P 2 '-5 s2 v 20 p j+ 4 p 2 1 

-[K s (L T -L T0 )+ c s i T ) l 5 If B H )/l t + Q *P 


(115) 


— O 

— (O 
r- rt 


— O 
mu 


11. 0 Equation 

E 

e p {I ypb } 


* 2 {(I xpb- I ypb ) ^“p 1 }+ Ira P r m ps ] { V p [ 5 S2 V 20 p ]+ Vp IC V 29 p ] 


+y p z p [ ’ s V"p— p 


c2e P ]+i p 2l "^ c2 v 2: p 1+ V I “2 s2 v 20 p ]+ V t ^ S2e p J} 


A 2 r l. 


-I K s (L T - L T0 )+c s f, T 1 rA H + B H n + c ^ /l t + Q ep 


(116) 
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SECTION 8 - SIMPLIFIED EQUATIONS Or MOTION 

Equations (106) to (116) are for the most general situation 
possible, and as a result, are quite lengthy. Under some 
circumstances, these equations can be simplified. If this can 
be accomplished, a significant decrease in computer time will 
be realized. The first simplification occurs if the forebody's 
aerodynamic (body) reference axes are principal axes. In this 
case X = Y = 1 = 0 and l xyb = l xzt) = I yzb = Equations (106), 
(107), (108), (112), (113), and (114) become: 


X{m) = [K s (L t -L tq )+C s L t ] [A/L t J+F x 

(117) 

Y(m} = [K s (L t -L to )+C s L t ] [B/L T ]+F y 

(118) 

Z(m} = [K s (L T -L T0 )+C a L T ] [C/L T ]-mg + F z 

(119) 

(s 2 e)-t [i yb (s 2 4»)-*-i zt> (c 2 4 >) 3 (c 2 e)}+0{(i yb -i zb ) 

<jces 2 *) ) + *{-i xb c 13 


* ^Ul xb -I yb (S 2 <}>)-I zb (C 2 ^)) (-S20)}+^{(I yb -I zb ) (-C 2 0S2<j))} 
+§$U I xb -(Iy b - I zb ) (C2<j>) ] (C9) }+0 2 { (I yb “I zb ) (^9S2<j>)} 

- tK s ( V L TO )+C sV C * If + 1 If ]/L T + % 


( 120 ) 


^yb^zb* (^9S20)}+0{I yb (C 2 4))+I zb (S 2 <{»)} 
^ Ul xb +(I yb" I zb ) (C2<p)] (-C0)}+0^{(I yb -I zb ) (S24.) } 

5 ' 2 < II xb- I yb (s2 ^ ) - I zb (c2<f)1( I S29)} 


- 3 >- 
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-I K S (L T- L TO )+C S i T 1 < A H + B ft +C ^/ L T + Q 0 (121) 

, ^- I xb C 13> + ^ I xb } 

- 4'0{ '[i X b+( I yb - I Z b ) CC24.) 3 (C6) }+^ 2 { (I yb -I zb ) (^C 2 es24») 
+§ 2 {(I yb“ I 2 b )( “I S2 ^ )} 

" rK s (L T" L TO )+C s^T J tA if + B if + ^ If 3 / 11 ! + Q (f> (122) 


The above six equations of motion have not only been shortened, 
but they also have been uncoupled in the translational 
accelerations making them easier to solve. The second 
simplification involves the decelerator degrees of freedom. 

If the added masses of the decelerator are ignored (m pJl = ro p S 
= nip) , Equations (109), (110), (111), (115), and (116) become: 


VV - [k s ( V l to )+c A ] (a/l t' 

Yp{m p }= - [Kg (L t “L to ) +C s L t ) [B/L t ] 
Z p {m p }= -[K g (L T -L T0 ) + C g L T ] [5/1*] 


+ F 

(123) 

xp 

+ F m 

(124) 

yp 

- m g + F 

(125) 

P ZP 
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*p {I xpb (s 26 p> +I ypb (c 2 e p>> - Vp^-W^V 1 

• [k s (l t' l to )+C 3 I 't 1 t5 H“ + 5 i/l t + Q 4,p 

P P 


(126) 


VW - *p 2{(I xpb- I ypb> <2® 2 V> 


-'WW+Vy 1 !A W- + B If + C ]/L T + Q 0p 

P P P 


(127) 


Like the simplified equation for the forebody, the decelerator 
equations have also shortened. Furthermore, they have 
completely uncoupled in the second derivatives making numerical 
integration easy. 
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SECTION 9 - GENERALIZED FORCES - AERODYNAMICS 

The nonconservative forces acting on the forebody are due to 
aerodynamics o The aerodynamics and the convention used in this 
report apply to che Space Shuttle Solid Rocket Booster (S.R.B.). 

If a different body is to be simulated, the aerodynamic coefficients 
and possibly the convention used to define them, would change. 

For the S.R.B., the aerodynamics are a function of roll angle, 
angle-of-attack, and Mach number. The angle-of-attack is 
measured from the total velocity vector to the positive longitudinal 
axis (X^) as shown in Figure 3. 

, (Y ' * + Z, ; 

a = tan X I ] (128) 

X b 

The normal force coefficient, C N , is in the plane formed by 
the velocity vector and the longitudinal axis, and is perpendicular 
to the longitudinal axis (X^) . The roll angle, <)>^, is then 
measured from the normal force coefficient to the Z^body axis. 

The axial force coefficient is defined as usual, positive in 
the negative X^ direction. Finally, the side force coefficient 
is perpendicular to the X^ body axis and to the normal force, 
such that the directions of C^, C N , Cy form a right-handed 
orthogonal coordinate system. Mathematically, the aerodynamics 
roll angle is given by: 

<J>^ = tan -1 [-Y b /-Z b ] (129) 


- 38 - 
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The positive directions of the moment coefficients are shown 
in Figure 3 as double arrows. Damping ant coefficients 
are about the body axes (X^, Y^, Z^) . Aerodynamic body axes 


forces are 

given as: 


F , = 

xb 

-q S C A 

(130) 

F . = 

yb 

qSl-CyC^ + C N S0 i 3 

(131) 

F , = 

zb 

qS[C y S$ i + C N C$ i J 

(132) 


The body axes forces are converted to inertial a : is force 
using the elements of [CJ , Equation 4. 


F x F xb C 11 + F yb C 21 + F zb C 31 


F y = F xb C 12 + F yb C 22 + F zb C 32 


F_ = F , + F. , + F , C 


xb C 13 + yb C 23 + F zb 33 “ mg 


Body axes torques are: 


(133) 

(134) 

(135) 


xb 


yb 


zb 


* gSd ^ + C £p ^ 


w xb d 

(- 58 --)] 


to , d 
vb 


■ ? Sd '- C m C *i - C n S *i + C m HSv"> ] 

q 


“zb d 


= qSdtC S*. - C n C*. + C n <-*$-» 


r 

The body axis torques are transformed to generalized 
torques using Equations (5) and (1) and noting sign 
conventions. 


(136) 

(137) 

(138) 


- 1 + 0 - 
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‘ tT xb C 13 + T yb°23 + T zb = 33 ! 

(139) 

« 

Q e a 

- T yb C4 + T *b S * 

(140) 


O 

-e- 

1) 

T xb 

(141) 



i 


The aerodynamics of a decelerator (parachute) are not well 
known because a parachute is not a rigid body, ana does not 
lend itself to easily obtainable test data, especially 
under dynamic conditions. Consequently the aerodynamics 
of a symmetric decelerator tend tc be relatively simple 
due to a lack of better understanding rather than the 
inability to use available information. If better 
aerodynamic data is attainable, it is a simple matter to 
alter the body forces and torques appropriately. 

For this report, the decelerator body forces and torques are: 


xpb 

— 

_q p S p C Ap 



(142) 

l 

ypb 

= 

q p S p C NP S< "pi 



(143) 

t 

zpb 


q S C x Cd> . 

^p p Kp v pi 



(144) 

1 

ypb 

= 

— q S d (C _ - 

^p p p V MP 

• 

o.i * e * d /v ) 

p p p 

c0 . 

pi 

(145) 

l 

zpb 


^p S p d p <C MP • 

• 

0.1 * 0 * D /V ) 

P P P 

S0 . 
Pi 

(146) 

Pi 

a 

tan-h-Yp^-^] 



(147) 
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The generalized forces for the decelerator are: 


xp 

F xpb 

c pll + 

F , 
ypb 

C p21 + 

F 

r zpb 

C p31 

1 — 
yp 

F xpb 

C pl2 + 

F , 
ypb 

C p22 + 

F . 
zpb 

C p32 

1 _ 

zp 

F xpb 

C pl3 + 

F . 
zpb 

C p33 

m g 
P 


1 . 

■'J'p 

T zpb 

C p33 





! 0p = 

-T , 

ypb 







(148) 

(149) 

(150) 

(151) 

(152) 
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CHAPTER III COMPUTER PROGRAM 

SECTION 1 - FEATURES OF THE COMPILER PROGRAM 

The computer program contains the following features. 

1. The program has many options which simplify the input of data or decrease the 

program run time. Use of the options are contained in the listing of the 

program as comment cards. These options are: 

a. Opti ms are included which change the dimensions of the aerodynamic 
coefficient arrays as dictated by input requirements. 

b. An option is provided (OPDT = 1.0) which automatically determines the 
magnitude of the integration time interval, DT. 

c. An option is provided (OPSP = 1.0) which calculates the parachute drag 
area (SP) time history. If this option is rot used the drag area versus 
time is input into the program in the form of look-up arrays. 

d. An option for including longitudinal and lateral added air mass effects 
on the parachute (OpAM = 1.0) is included in the program. 

e. A provision is made to use simplified equations of motion (0P0S = 0.0) to 
reduce run time, if all the forebody products of inertia and center of 
mass offsets are equal to zero. 

f. An opti n (OPPLCfT = 1.0) for making a plot tape is available. 

g. English or metric systems may be used for data input and out by equating 
0METRC to 0.0 or 1.0 respectively. 

2. All aerodynamic coefficients are read into the program as functions of angle 

of attack, roll angle, and mach number in the form of three dimensional look-up 

arrays . 
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3. The initial start conditions for the forebody and aft body are completely • 

general. 

4. The stacking of design cases is possible. 

5. The attachment location of the tether to the forebody is completely general. 

6. The tether load and the angle it makes with the centerline of the forebody are 
program outputs. 

7. All load and trajectory data are output at pre-selected times. 

8. Termination of a design case occurs at a predetermined time or altitude. 

9. The program calculates the effective system spring constant. 

10. The program calculates the parachute physical properties as the parachute 
inflates as a function of time. 

11. The parachute may have three stages of reefing, if the automatic drag area 
versus time option is chosen. 

12. As the parachute inflates, the drag area versus time follows a second degree 
curve (y = ax^) . 
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SECTION 2 - INPUT 

Except for the variable COM, all inputs are read in under the format statement 
8F10.0. COM is an 80 column header card. All of the following variables are 


defined in 

the nomenclature. 




INPUT ITEM 

VARIABLE 

NUMBER OF CARDS 

a) 

AIPHI , AIPHID, AJALPF, AJALPM, ARAM, AKAMD, 





OPSYM, OPDA 

1 

card 


b) 

PPHIE 

1 

card 


c) 

AALPFE 

1 

or 2 

cards 

d) 

AALPME 

1 

or 2 

cards 

e) 

AAM 

1 

card 


f) 

CCA 

4 

to 128 cards 

g) 

CCN 

4 

to 128 cards 

h) 

CCLM 

4 

to 128 cards 

i) 

CCY 

0 

or 4 

to 128 cards 

j) 

CCLL 

0 

or 4 

to 128 cards 

k) 

CCLN 

0 

or 4 

to 128 cards 

1) 

CLIP, CLMQ, CLNR 

0 

or 1 

card 

m) 

PPHIDE 

0 

or 1 

card 

n) 

AALPDE 

0 

or 1 

card 

o) 

AAMD 

0 

or 1 

card 

P> 

CCLLP 

0 

or 4 

to 64 cards 

q) 

CCLMQ 

0 

or 4 

to 64 cards 

r) 

CCLNR 

0 

or 4 

to 64 cards 

s) 

AALPPE 

1 

card 


t) 

AAMP 

1 

card 


u) 

CCAP 

2 

to 8 

cards 

v) 

CCNP 

2 

to 8 

cards 

w) 

CCMP 

2 

to 8 

cards 
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INPUT ITEM VARIABLE NUMBER OF CARDS 


x) 

PS, PT, EPL, EPT 

1 

card 

y) 

X, Y, X, V, GAME, CHIE, EPSI, ETA I 

1 

card 

z) 

WT, IXB, IYB, IZB, IXYB, IXZB, IYZB 

1 

card 

aa) 

S, D, XBAR, YBAR, ZBAR, OPPRIN, OPPLOT, OPOS 

1 

card 

bb) 

PSIE, THEE, PH IE , OMXBE , OMYBE, OMZBE 

1 

card 

cc) 

A, B, C, OPAM, OMETRC 

1 

card 

dd) 

PSIPE, THEPE, PSIPDE , THEPDE , VP, GAMPE , CHIPE 

1 

card 

ee) 

LS, LTO, DLTO, NS, NT, DP, CCRIT 

1 

card 

ff) 

AMAX1 , AMAX2 , DSX1, DSX2 , AMAY1 , AMAY2 , DSYl, DSY2 

1 

card 

gg) 

WTC, WTL, OPSP, OPDT 




IF (OPSP .EQ. 0.0) GO TO ITEM kk) 



hh) 

TRO, TR1 , TR2 , TR3 

1 

card 

jj) 

SPRO, SPRl , SPR2 , SPR3 , PCT01, PCT02 , PCT03 , POROS 

1 

card 

kk) 

TTIP , SSP 

4 

cards 


IF (0PDT.EQ.0.0) GO TO ITEM mm) 



ID 

GLOAD, FSULT, AERATO, TO, DTP1, TTT, HHH 

1 

card 

mm) 

DTI, TO, DTP1 , TTT, HHH 

1 

card 

nn) 

COM 

1 

card 


The values of the variables read in input item "a" determine, in part, the sizes 
of the aerodynamic arrays. The axial and moment coefficients have the option of 
using either eight or sixteen angles-of-attack (one or two cards). If, for example, 
five or eleven angles-of-attack are needed, one or two cards are needed respec- 
tively. The roll and Mach number arrays may vary from two to eight. As an 
example consider the array CCA where the value of C^ depends on five roll angles 
C^), eleven angles-of-attack (0( ) and seven Mach numbers (AM). The array 
PPHIE would be read in on one card containing five distinct roll angles, the last 
three fields of ten digits would be blanks. The array AAU>FE would be read in on 

two cards. The first card would contain eight distinct angles-of-attack, and the 

second card would contain three distinct angles-of-attack and five blank fields cf 
ten digits. The array AAM would be read in on one card containing seven distinct 

Mach numbers and one blank field of ten digits. The first element in each of the 

above arrays should start at zero and increase numerically until ♦ c highest 
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possible value expected to be encountered is specified. In this particular example, 
the array size used will be CCA(5, 16, 7). The proper read sequence is to first 
read two cards containing the values of at eleven angle -of-attack, the initial 
roll angle (zero) and the initial Mach number (zero). These cards are followed by 
two cards containing values of at eleven angles-of-attack, the second roll angle 
and the initial Mach number. This is continued for five roll angles at the initial 
Mach number. After these ten cards, the same procedure is followed for the s :ond 
Mach number, and the third, etc. up to seven sets of ten cards. 

All the aerodynamic coefficient arrays are read similarly. However, notice that 
the angle-of-attack array associated with the moment coefficients is different than 
that associated with the force coefficients. Also, the damping moment coefficient 
arrays (input items "p"» "q", ana "r") may not be read in at all, depending on the 
value of OPDA. Instead, input item "1" can be used if the damping coefficients are 
constant. Finally, the damping coefficients correspond to the arrays read in input 
x terns m j n j 3 nd o * 

Figure 4 helps co clarify the meaning of the input parameters associated with the 
added air mass on the parachute (Ref. Input Item ff). 



Parachute Reference Diameter (ft) (log) 


Figure 4. Input Parameters for Parachute Added Air Mass 
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Figure 5 helps to clarify the meaning of the parameters associated with OPSO = 1.0 
which directs the program to calculate the parachute area time history as time 
ad 'ances (Reference Input Items hh and jj). 
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Figure 5. Parameters for Program Calculated Parachute Drag Area 
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SECTION 3 - OUTPUT 


All output variables are defined in the nomenclature. Before beginning the simu- 
lation, the following variables, specifying the characteristics of the rigid body 
and initial parameters, are printed out. 

Line 1. COM 

Line 2. IXB, IXYB, 1XBAR, S, CLIP, OPPRIN, OPSYM, AIPHI, AIPHID, DTI EPSI 

Line 3. IYB, IXZB, YBAR, D, CLMQ, 0PPL0T, OPDA, AJALPF , AJALPM, AJALPD , 

TTT, ETA1 

Line 4. IZB, IYZB, ZBAR, WT, CLNR, OPOS, OMETRC , ARAM, AKAMD, HHH 


If CLLP , CLMQ, and CLNR are constants for the simulation, their values are printed 
out in the appropriate place. If the damping coefficients are found from interpo- 
lation of three dimensional arrays, CLLP, CLMQ, and CLNR are set equal to zero for 
this printout only. Several variables dealing with the decelerator are then 
printed out. 


Line 5. 


Line 6. 


Line 7. 
Line 8. 


A, LTO, LS, AMAXl , AMAX2 , AMAYl , AMAY2 , AP, GLOAD, FREQP, OPAM, 
PCT01 

B, NT, NS, DSK1, DSK2 , DSY1, DSY2 , CHIPE, FSULT , POROS, OPDT, 
PCT02 

C, DLTO, DP, WTC, WTL, WTP, CCRIT, VP, AERATO, TO OPSO, PCT03 
TRO, TR1 , TR2 , TR3 , SPRO, SPR1, SPR2 , SPR3 


The parachute suspension line load and strain arrays are printed out next on Lines 9 
and 10, 


Line 9. PS(I) 

Line 10. EPL(I) 


The tether line load and strain arrays are printed out next on Lines 11 and 12. 

Line 11. PT(I) 

Line 12. EPT(I) 
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The parachute inflation time history array and drag area array are orinted out 
next. If the option (OPSP * 1.0) the arrays are set equal to zero Lc>:ause they 
are not known before initial time TO. 

Line 13. TTIP(I) 

Line V. SSP(I) 

The computer program then checks the option variable OPPRIN. If OPPRIN = 1. , 
all the aerodynamic data is listed as follows: 

PPHIE(I) 

AALPFZ (J) 

AALPME (J) 

AAM(K) 

CCA(I,J,K) 

CCN(I,J,K) 

CCLM ( I , J , K) 

If OPSYM = 0., the following aerodynamic data is listed 
CCY(I,J,K) 

CCLL(1,J,K) 

CCLN(I,J,K) 

In the above aerodynamic coefficient arrays, AALPFE(J) is associated with CCA, 
CCN, and CCY; AALPME(J) is associated with CCI^l, CCLL, and CCLN. 

If OPDA = 1., the damping aerodynamics is listed. 

PPHIDE (I) 

AALPDE (J) 

AAMD(K) 

CCLLP(I, J,K) 

CCLMQ(I,J,K) 

CCLNR(I, J,K) 


- 50 - 


"~\~7 


REF: EOI 380 


I 


GOODYEAR AEROSPACE 

COMPOM AXIOM 

GER- 16047 


The aerodynamic arrays associated with the decelerator then follow. 

AALPPE (1) 

AAMP(I) 

CCAP(l) 

CCNP(I) 

CCMP(I) 

After the listing of the input data, the computer program begins numerically inte- 
grating. At T = T q and at predetermined time increments, the following data is 
printed out. 

Line 1. T, X, XD, XDD, FX, CA, V, TENS, XP, XPD, XPDD, FXt, CDAP, CMP 

Line 2. TXB, Y, YD, YDD, FY, CN, AM, LT, YP, YPD, YPDD, FYP, CNP, AMP 

Line 3. TYB, Z, ZD, ZDD, FZ, CY, DYPR, TPD, ZP, ZPD, ZPDD, FZP, TYPB, DYPRP 

Line 4. TZB, PSIE, PSIDE, PSIDDE, QPSI, CLN, ALPE, OMXBE, PSIPE, PSIPDE, 

PSPDDE , QPSIP , TZPB, ALPPE 

Line 5. GAME. THEE, THEDE, THEDDE , QTHE, CLM, PHIIE, OMYBE, THEPE, 

THEPDE , THPDDE , QTHEP , TPDXB , GAMPE 

Line 6. CHIE, PHIE, PHIDE, PH^DDE , QPHI, CLL, PHIAE , OMZBE, KS, CLLP, 

CLMQ, CLNR, TYPDRB , PULAN 

Line 7. MPAL, MPAS, DMD, QMAXPB, 1XPB, IYPB, SPD, SP, SPRU, SPRL, TINT, 

TNINY, TFI , XPBDI 

When the simulation reaches HHH or TTT, the computer will write out "RUN ENDED 
BY CONSTRAINTS". It will then attempt to read in more data cards, to initialize 
for another simulation, starting with input item "y". If there are no data 
cards available, the program will CALL EXIT. 
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SECTION 4 - NUMERICAL SOLUTION 

For the most general type rigid body, there are six second order differential 
equations, coupled in the acceleration terms. These six equations can be written 
as : 

6 

Z D i;J U. = 3 = 1,2,..., 6 (153) 

i=l 

and solved simultaneously using the PIVERT subroutine. PIVERT uses Gauss elimina- 
tion with complete pivoting to obtain the largest diagonal elements. After solving 

00 

for the accelerations (U. in equation (153)) j the results are numerically inte- 

1 ( 3 ) 

grated using Runge-Kutta, fourth order techniques' . 

If the forebody has the properties that X = Y - Z = I ^ = I xz ^ = 1^.^ = 0. , the 

equations of motion greatly simplify for the forebody. In the case of integrating 
the Euler angles, three equations remain coupled in the acceleration terms, and are 
separated using PIVERT. The three translational accelerations are already in a 
suitable form to integrate immediately. A simpler situation occurs if the added 
masses of the decelerator are neglected. All five equations of motion are un- 
coupled in the second derivative and are easily integrated by lith order Runge- 
Kutta . 
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SECTION 5 - PLOTTING ROUTINE 

If 0PPL0T = 1., eleven variables are saved in arrays. At the end of the simula- 
tion, any or all of these variables are plotted by calling PLTRAJ and setting the 
appropriate arguments. PLTRAJ was originally written for use on a CALCOMP 563 
plotter and 750 tape drive. It has been modified for use at M.S.F.C. where a 
SC 4020 plotter is the preferred plotter. The original PLTRAJ will plot up to 
4 variables versus time on one graph for each call to PLTRAJ. The modified PLTRAJ 
for the SC 4020 plotter plots only one variable versus time per plot; therefore 
four plots will be made instead of one for each call to PLTRAJ. Two hundred data 
points are plotted on each graph per variable. 
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SECTION 6 - ENGLISH TO METRIC CONVERSION 

The computer program operates in either English or Metric units, '’’he program input 
and output is in English units unless the option parameter, OMETRC, is set equal to 
1. If OMETRC = 1. the input and output is in the Metric System. A conversion 
table from English to Metric is given belov for commonly used engineering parameters. 


ENGLISH TO METRIC CONVERSION 
REFERENCE NASA SP 7012 
♦EXACT 


FJRCE 

(LB) X 4.4482216152605* = (1) NEWTON 

N 

LENGTH 

(FT) X .30480000 * = (1) METER 

m 

MASS 

(SLUG) X 14.5939029 = (1) KILOGRAM 

kg 

SPEED 

(FT/SEC) X .3048 = (1) METERS SEC 

m/sec 

PRESSURE 

(lb/ft 2 ) x 47.880258 = (1 ) newton/meter 2 

N/m 2 

volume 

(FT3) x . 028316846592* = METERS 3 

m 3 

AREA 

(FT 2 ) X .09290304* = METERS 2 

2 

m 

ACCELERATION 

(ft/sec 2 ) X .3048* = meter/sec 2 

m/sec' 

INERTIA 

(SLUG-FT 2 ) X 1.355817945-- KILOGRAM-METER 2 

, 2 

kg-m 

TORQUE 

(FT -LB) X 1.355817948 = METER - NEWTON 

m-N 

DENSITY 

(slug/ft 3 ) X 515.379 = kilogram/meter 3 

kg/m 3 

viscosity 

(SLUG/FT-3EC) X 47.880258 = NEWTON SEC/METER 2 

(N-sec'/i 

SPRING CONSTANT 

(lb/ft) X 14.59390293 = newton/meter 

N/m 


* Exact Numbers - No round offs 


GOODYEAR AEROSPACE 

CORPORATION 

GER-1^047 


SECTION 7 - CONCLUSIONS AND RECOMMENDATIONS 

The following conclusions and recommendations are made for use of the 6+3 D.O.F. 
computer program. 

1. The (6+5) DOF loads assessment computer program should be used primarily to 
analyze the loads induced on a wobbling or spinning body when the body is 
stabilized by the deployment of a drogue parachute. 

2. After the body has been stabilized by the drogue, further parachute deployments 
(main chutes) should be analyzed using th planer (3+3) DOF computer program. 

The (3+3) program should be used because of the following reasons: 

a. The (3+3) program is faster and easier to use than the (6+5) program. 

O 

b. Terminal descent with .he forebody pitch angle equal to ±90. presents 
no mathematical solution problem using the (3+3) DOF program. 

3. It should be noted here that the (6+5) DOF program has a mathematical singularity 
point at a forebody pitch angle of ±90°. To permit passage through this point 

the six forebody accelerations are frozen at their last value when the pitch 
angle is in the region of 89,8° < " < 90.2°, This ' iduces some error in the 
translational coordinates and attitude of the forebody, but it has been shown 
to be small for normal velocity passes througn this point. A time count.e 
(TNINY) for the time spent in this region is a program ovtput . 
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CHAPTER IV - PROGRAM LISTING AND SAMPLE COMPUTER RUN 

The following program listing is for the Uni vac 1108 at M.S.F.C. and adapted from 
che IBM 360 listing used by Goodyear Aerospace Corporation. 

The sample problem, SRB Stabilization by 54' Drogue Paiachute, represents the 
deployment of a 54' drogue from a SRB which is wobbling and flying broadside 
to the wina vector. The trajectory of the SRB is nearly vertical. The drogue 
is starting to inflate and is stretched out normal to the SRB centerline. The 
drogue has one stage of reefing (0.82 of full open area). 

Some of the more important initial conditions are given in the table below. 


Atltitude 

Ft 

19,000. 

Velocity 

Ft/Sec 

553. 

Angle of Attack Deg. 

90. 

Flight Path Angle Deg. 

85. 

Body Axis 

Rates 


Pitch 

Deg/Sec 

1.0 

Yaw 

Deg/Sec 

0.0 

Roll 

Deg/Sec 

0.0 


The output from the sample problem starts on Page 57, and selected portions of 
the trajectory are found starting on page 197. 
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